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Abstract 

The class of the two-loop massless crossed boxes, with light-like external legs, is the 
final unresolved issue in the program of computing the scattering amplitudes of 2 — > 2 
massless particles at next-to-next-to- leading order. 

In this paper, we describe an algorithm for the tensor reduction of such diagrams. 
After connecting tensor integrals to scalar ones with arbitrary powers of propagators 
in higher dimensions, we derive recurrence relations from integration- by-parts and 
Lorentz-invariance identities, that allow us to write the scalar integrals as a combination 
of two master crossed boxes plus simpler-topology diagrams. 

We derive the system of differential equations that the two master integrals satisfy 
using two different methods, and we use one of these equations to express the second 
master integral as a function of the first one, already known in the literature. We then 
give the analytic expansion of the second master integral as a function of e = (4 — D) /2, 
where D is the space-time dimension, up to order O (e ) . 



■"■e-mail: Ch. AnastasiouOdurham. ac .uk 

2 e-mail: Thomas . GehrmannOphys . uni-karlsruhe . de 

3 e-mail: Carlo . Olear i@durham .ac.uk 

4 e-mail: Ettore . RemiddiObo . inf n. it 

5 e-mail: tauskOphysik . uni-f reiburg . de 



1 Introduction 



The increasing precision of high-energy scattering experiments has reached the level that 
next-to-next-to-leading theoretical results for 2 — >■ 2 scattering processes of massless particles 
are demanded |T[. Several steps have already been taken towards this goal: in Refs. |3| 
the tensor reduction of two-loop planar boxes and the computation of the relevant master 
integrals have been given. Reference Q dealt with the tensor reduction of the pentabox, 
while other simpler two-loop diagrams obtained by one-loop insertion into one-loop four- 
point functions have been treated in Ref. ||. 

The final unresolved issue is the tensor reduction of two-loop crossed boxes of Fig. [I] and 
the evaluation of the corresponding master integrals. One of the master integrals, corre- 
sponding to the scalar integral with all powers of propagators equal to unity, was computed 
in Ref. || as an analytic expansion in e = (4 — D)/2, where D is the space-time dimension. 

A first result of this paper is to present a method to reduce tensor integrals of crossed two- 
loop boxes to two crossed-box master integrals, plus master integrals for simpler-topology 
diagrams. We choose as first master integral the already evaluated crossed box with all 
powers of propagators equal to one, and as second master integral the scalar integral where 
the power of the second propagator is equal to two, all the others being one. 

As it is well known j|, [7], tensor integrals can be related to scalar integrals with higher 
powers of the propagators in higher dimensions. We make use of recurrence relations obtained 
by integration-by-parts || and Lorentz-invariance identities to connect integrals with 
different powers of the propagators and to derive a reduction algorithm that expresses the 
generic scalar integral as a function of the two master crossed boxes, plus simpler sub- 
topologies. We explicitly give the equations that connect the two master integrals in different 
dimensions (dimensional shift). 

Quite in general, it has been shown that the master integrals for any Feynman graph 
topology satisfy a system of first-order differential equations on any of the Mandelstam 



variables on which they depend ||10| . We derive a system of two coupled differential equations 
for the crossed-box master integrals by two independent methods: first using the raising and 
lowering operators and second by taking the on-shell limit of the differential equations for 
the crossed box with one off-shell leg ||. 

Inserting the analytic expression of the first master crossed box into the differential 
equation for the first master integral, we obtain an algebraic equation for the second master 
crossed box, that can be solved to give the analytic expansion in e of the second master 
integral. The differential equation for the second master integral provides an overall check 
of the calculation. 
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Our paper is organized as follows: in Section |2] we introduce the notation we are going 
to use. In Section |3] we briefly review how tensor integrals are related to scalar ones, and 
we give the reduction algorithm to express these integrals as a combination of the two 
master crossed boxes, plus simpler master integrals for topologies with fewer propagators. 
In Sections ^ and |5] we derive the differential equations that the two master integrals satisfy 
in arbitrary dimensions D, using the two different methods. The analytic expansion of the 
second master integral is given in Section || For completeness, in Appendix |A] we present the 
algorithm for the tensor reduction and for the dimensional shift of the crossed triangle, since 
this is one of the non-trivial sub-topologies produced by the reduction procedure. Finally 
we conclude with Section [7]. 




We denote the generic two-loop tensor crossed (or non-planar) four-point function in D 
dimensions of Fig. [I] with seven propagators Ai raised to arbitrary powers Ui as 



Xbox D (vx, u 2 , i/ 3 , 1/4, us, u 6 , u 7 ; s, t) [1; P; P; Wk"; fcT; . . .] 

d D k r d D l [1; P; k»k v ; kH h 



where the propagators are 



i7r D/2 J m D/2 A^AfA^A^A^A^A? ' 



A x = (k + I + p u f + iO, 
A 2 = (k + l+p 1M ) 2 +i0, 
A 3 = (k + l) 2 + i0, 



(2.1) 



3 



A 4 = l 2 + i0, (2.2) 

A = (l+p 3 ) 2 + i0, 

Aq = k 2 + iO, 

A 7 = (k + p 4 ) 2 + iO. 

The external momenta pj are in-going and light-like, pj — 0, j — 1 . . . 4, so that the only 
momentum scales are the usual Mandelstam variables s = (pi + P2) 2 an d £ = Q»2 + P3) 2 , 
together with u — — s — t. For ease of notation, we define Pij = Pi +pj and p^ = Pi +pj +Pk- 
In the square brackets we keep trace of the tensor structure that may be present in the 
numerator of Eq. (|2.1| ) . When the numerator is unity we have the scalar integral 

Xbox D (u u u 2 , v 3l ^4, ^5, "6, vf, s, t) [1] = Xbox D (u u u 2 , u 3 , 1/4, z/ 5 , u 6 , u 7 ; s, t) . (2.3) 



3 Tensor reduction using raising and lowering operators 

Tensor integrals can be related to combinations of scalar integrals with higher powers of 
propagators and/or different values of D J7|, f§. In fact, introducing the Schwinger parameters 
Xi, we can write the integrand of Eq. (|2.1|) in the form 



where 



and 



A ^ 1 A u 7 = J Vx ex P (j2 x i A ^j , (3- 1 ) 



Vx = U { -FTVl dx iX r\ (3.2) 

i=l 1 \ u i. 







^XiAi = ak 2 + bl 2 + 2ck-l + 2d-k + 2e-l + f, (3.3) 



i=i 



a = xi + x 2 + x 3 + xq + x 7 

b = xx + x 2 + x 3 + x 5 + x 4 

C = Xi + x 2 + x 3 

# = ^34 + ^134 + ^4 

e M = ^34 + ^134 + ^3 

/ = x lS . (3.4) 
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We can diagonalize the exponent with the change of variables 



cL \ M 

fc" -> [K + X) , (3.5) 



a 

F -> (L + 3T, (3.6) 



where 



X^ = ^| - [(x 2 + xi)(x 7 + x 5 + x 4 ) + (x 5 + x 4 + x 3 )x 7 ]p 4 l + [x 3 x 5 - x 4 (x 2 + Xl)]P3 

- x 2 (x 5 + x 4 )p^, (3.7) 

y = ^{ [x 3 x 7 - x 6 {x 2 + xi)]p% - [(x 2 + x x )(x 7 + x 6 + x 5 ) + (x 7 + x 6 + x 3 )x 5 )]p^ 

-x 2 (a; 7 + a; 6 )rf}, (3.8) 



and 



P = (rr 7 + x 6 + a; 5 + x 4 ) (x 3 + 0:2 + £1) + (x 5 + x 4 ) (x 7 + x 6 ) . (3.9) 



The generic tensor integral (dropping all the dependences on vi and on the external scales) 
becomes 

d D K r d D L 



Xbox^P 1 . . . k^r . . . r n ] = jvx J J 



in D/2 

<-L \ w /„ cL 



k - ^ + x) ...(k-^ + x) (l + yr + yr 

(a^ 2 + ^ 2 + |), (3.10) 



x cxp 
where 

Q = x 2 (x 5 x 6 - x 4 x 7 )t + [x 1 x 3 (x 7 + Xq + x 5 + x 4 ) + x 3 x 5 x 7 - x 2 x 4 x 7 + XiX 4 X b ] s. (3.11) 

The integration over the loop momenta K and L is now straightforward. For example, the 
scalar integral and the two one-index tensor integrals are given by 

Xbox D is 2 , v 3 , v 4 , z/ 5 , u 6 , v 7 \ s,t) = Jvxl, (3.12) 

Xbox D (1/1, v 2 , v 3 , v 4 , u 5 , z/ 6 , v 7 \ s, t) [A; M ] = J Vx X^X, (3.13) 

Xbox D (is u i/ 2 , u 3 , v 4) z/ 5 , u 6 , vi\ s, t) [F] = jvx yi, (3.14) 

where 

expf^V (3.15) 



y<p 
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Recalling the definition ( |3.2| ), we see that we can absorb the factors X{ of X^ and into 
T>x, increasing the power of the i-th propagator by one 

1 i ( I * =► -v ± \> ** ee (3.16) 

where 

^Xbox (...,!/<,...) = Xbox D (. . . , i/j ± 1, . . .) , (3.17) 
while the factor V can be absorbed into X (see Eq. Q3.15 )) 
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(3.18) 



J>D/2 <p <p(D+2)/2> 

so that 1/V acts dimension increaser 

I =► d + , P =► d" (3.19) 

where 

d ± Xbox D = Xbox D±2 . (3.20) 
For example, using the expression for X^ of Eq. (|3.7| ), we can write 

Xbox D [P] = { [z/33 + z/55 + -z/44 + (z/ 2 2+ + z/ 1 l+)]^-z/ 2 2+(z/ 5 5 + + z/ 4 4+)^ 

[u 2 2 + + vil + ) (v 7 7 + + z/ 5 5 + + u 4 4+) 

+ (z/ 5 5+ + z/ 4 4 + + u 3 3+) v 7 7 + ]p>Z} Xbox D+2 [l]. (3.21) 

The task to compute tensor integrals has then been moved to the computation of scalar 
integrals with higher powers of the propagators in higher dimensions. 

In the next sections, we will follow the procedure already used in Ref. ||: 

1) we first reduce the powers of the propagators of the generic scalar diagram and we 
express it in terms of a finite set (basis) of master diagrams; 

2) then we relate the master integrals in higher dimensions to the master integrals in D 
dimensions (dimensional shift). 



3.1 The scalar crossed-box reduction 

The strategy to reduce the generic scalar integral to a linear combination of master ones is 
based on identities that relate scalar integrals with different powers of propagators. Some of 
these identities can be obtained using the integration- by-parts method ||. Since there are 
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three external independent momenta and two loop momenta, we can build ten identities, 
imposing that 

r rl D h r rl D l f) 

/^/^^l"/(M, P ,)]=», (3-22) 



where 



of = k", P (3.23) 

6" = V 1 , P, rf, (3.24) 
1 

AVi Al>2 A^3 AV4 4^5 A^B A u 7 

71 1 ^2 ^3 ^4 ^5 ^6 ^7 



f{k,l,Pi) AU 2 AU 3 AU 4 A U 5 A U 6 AU 7 ■ (3.25) 



Not all the scalar products that appear in the application of the integration-by-parts iden- 
tities can be written in terms of combinations of propagators: with the two loop momenta 
and with the three linearly independent external ones, we can form nine scalar products 
involving k and I. Since we have seven linearly independent propagators, we are left with 
two irreducible scalar products in the numerator, that we choose to be (I • pi) and (I ■ pi). 

Three more identities among the scalar integrals can be derived if we exploit the Lorentz 
invariance of the Feynman diagrams |J. In fact, since the Feynman integral is a function 
only of scalar products of the external momenta, it is invariant under the (infinitesimal) 
rotation 

p't = A^ u Pi, A Mi , = g^u + S e^, = -e Ufl . (3.26) 

We can then write 

dPk r d D l P/I , , r d D k r d D l r/1 , .. . „. 

f{k,l,pi) = / -™ / —f^ mi^. (3.27) 



m D/2 J m D/2 JVT'iriJ j mD/2 j 

Expanding in a Taylor series in 5 the right-hand side of Eq. (|3.27| ), we obtain 

d D k r d D l * df(k,l,pi) „ „ * \ 

Using the three independent external momenta, we can build three independent antisym- 
metric tensors 





U V 

= P1P2 


U V 

-P2P1 


€ 2 


= P1P3 


-P3P1 


av 
4 


u v 
= P2P3 


u v 
-P3P2 



that, once inserted into Eq. (|3.28 ), give rise to three more identities. 



Taking linear combinations of these thirteen equations, we can build the following system 
sv x \+ - z/ 7 7+6 _ - z/ 5 5 + 4~ - (z/ 2 2+ + ^1+) 3 _ - u 1257 - 2z/ 346 + 2D = (3.29) 

7 



sz/ 3 3+ - z/ 6 6+7- - z/ 4 4+5- - (V 3 3+ + z/ 2 2+) 1" - z/ 2346 - 2z/ 157 + 2D = (3.30) 

2 (Z • p 4 ) ^4+ - (z/ 6 6 + + z/ 5 5 + ) 7 _ + z/ 5 5 + l~ + z/ 4 4 + (3 _ - 6 _ ) 

- z/ 456 - 2z/ 7 + D = (3.31) 
2 (Z • p 4 ) + z/ 7 7+6 _ + z/ 5 5+ (7- - l - + s) + z/ 4 4+ (6 _ - 3 _ ) 

+ ^457 + 2z/ 6 - D = (3.32) 
2 (Z • p 4 ) ^6+ + z/ 7 7+ (5- - l - ) + z/ 6 6+ (7- + 5 _ - l - + s) + z/ 4 4+5 - 

+ z/ 47 + 2z/ 5 - £> = (3.33) 
2 (Z • p 4 ) ^7+ + z/ 6 6 + (3- - 4 _ ) - z/ 7 7+ (6 _ + 4- - 3 _ ) - v 5 5+4~ 

- z/ 56 - 2z/ 4 + £> = (3.34) 
2s (Z • p 4 ) u 2 2 + + 4s (Z • p 4 ) z/ 3 3+ + z/ 2 2+ [(* + s) (s + 1") - t3~ 

+ s (2z/ 6 6+ + 2z/ 3 3+ + z/ 2 2+) 7" - s (2i/ 3 3+ + z/ 2 2+) 6" 

- (2D - 2z/ 13457 - v 2 ) s = (3.35) 
2s (Z • p 4 ) z/ 3 3+ - 2s (Z • p 4 ) ^1+ + (t + s)z/ 2 2+ (l~ - 3~) + s (z/il + - i/ 3 3 + ) 6~ 

+ s (z/ 6 6+ + i/ 3 3+ - z/il + ) 7 _ - sz/ 5 5+4 _ - sz/xl+3 - 

- (^6 - ^15) s = (3.36) 
2 (Z • p 4 ) ^3+ - 2 (Z • pi) z/ 2 2+ + (t + s) z/ 2 2+ + (z/ 6 6+ + u 3 3+ + z/ 2 2+ + ^1+) 7 - 

- (z/ 2 2+ + z/xl + ) 5- - z/ 3 3+4- + z/ 1236 + 2z/ 7 - D = (3.37) 

2s (Z • p 4 ) z/ 3 3+ - 2s (Z • pi) z/ 7 7+ + z/ 7 7+ [(* + s) (s - 1~) + t5~ + s2~ 

+ z/ 6 6+ (t4T - t3~ + s7~) + sz/ 3 3+ (7~ - 6~) + (t + s)z/ 2 2+l~ + sz/il+2 - 

+ s (z/ 126 + 2z/ 3457 - 2D)-t(D- v 67 + u 2 - 2v 45 ) = (3.38) 
2s (Z • p^ z/ 3 3+ - 2s (Z • Pl ) z/ 5 5+ + (t + s) z/ 5 5+ (7 _ - 1~ + s) + (t + s)z/ 2 2+l- 

+ (t + s)z/ 4 4+ (6~ - 3~) + s (z/ 6 6+ + z/ 3 3+) 7~ - sz/ 3 3+6~ - sz/il+2 - 

+ s (^1456 - ^2 + 2z/ 7 - £>) + t (2z/ 67 + z/ 45 - i/ 2 - D) = (3.39) 
2s (Z • p 4 ) z/ 3 3+ + 2s (Z • pi) z/ 6 6+ + z/ 6 6+ [t (4~ - 3~) - s (2~ - 1~ + s)] - sz/ 4 4+5~ 

+ sz/ 3 3+ (7- - 6 _ - 2 - ) + tu 7 7 + (5- - l - ) + (t + s)u 2 2+l~ 

+ s (z/ 34 - i/ 2 ) + t {y&i + 2z/ 45 - i/ 2 - D) = (3.40) 
2s (Z • p 4 ) z/ 3 3+ + 2s (Z • p^ z/ 4 4+ + z/ 4 4+ [(* + s) (6~ - 3~) + s5~] + 2sz/ 6 6+7~ 



+ (t + s)v 5 5+ (7- - 1") + (t + s)u 2 2+l- + su 3 3+ (7- - 6- + 2") 

+ s (z/ 23 + 2z/i46 + 3z/ 5 + 4z/ 7 - 3D) + t (i/ 45 — i^a + 2z/ 67 - D) = 0, (3.41) 

where we have introduced the shorthand z/jj = + z/j, fyk = v% + z/j + z/&, etc. Each equation 
acts on the integrand of the generic crossed box before any loop integration has taken place. 

Equations ( |3.29|) and (|3.30| ) of the system, being independent of the two irreducible scalar 
products, need no further manipulation, and can be rewritten in the form 



8Vl l + = u 7 7+6~ + u 5 5+4~ + [u 2 2+ + u 1 l+)3- + u 1257 + 2u 346 -2D, (3.42) 
sz/ 3 3+ = z/ 6 6 + 7~ + z/ 4 4+5~ + (z/ 3 3+ + z/ 2 2+) 1" + z/ 2346 + 2z/ 157 - 2D. (3.43) 



By repeated application of these two identities, we can reduce V\ and z/ 3 to one. During this 
process, the generic scalar box is expressed as a linear combination of crossed-box diagrams 
with V\ = z/ 3 = 1 and diagrams belonging to simpler topologies, that originate when powers of 
propagators are reduced (pinched) to zero by the decreasing operators. We will deal with the 
pinched diagrams later, concentrating now on the reduction of the remaining propagators. 

In order to use the other equations of the system, we have to eliminate the irreducible 
scalar products in the numerator. 

For example, applying the operator z/ 7 7 + to Eq. Q3.31 ) and z/ 4 4 + to Eq. ( 3.34f) , and taking 
the difference, we get 



(D - 2v 7 - z/ 56 - 2) u 7 7 + — (D — u 56 - 2z/ 4 - 2) z/ 4 4+ - {y 7 - u A ) (z/ 5 5+ + z/ 6 6" 

+ z/ 5 5+z/ 7 7+l- - z/ 4 4+z/ 6 6+3- = 0. (3.44) 

In the same way, we can apply vqQ^ to Eq. ( |3.32| ) and z/ 5 5 + to Eq. ( |3.33| ) and take the 
difference, to obtain 



(D - z/ 47 - 2z/ 5 - 2) z/ 5 5+ — (D — v A7 - 2u 6 - 2) z/ 6 6+ + (u 6 - z/ 5 ) (z/ 4 4+ + u 7 7+ j 

+ i/ 5 5 + ^ 7 7 + l _ - i/ 4 4 + ^ 6 6 + 3 _ = 0. (3.45) 
Combining Eq. ( [3 .44 ) and ( |3.45|) to eliminate z/55 + , we have 

, 4 4 + = {D - 2 ^- 2 \ v 7 7+-2 f , 6 6+ 
(D — 2^45 — 2) (D — 2^45 — 2) 

z/ 5 5 + z/ 7 7 + l _ - z/ 4 4 + z/ 6 6 + 3 _ ) , (3.46) 



(D - z/ 4567 - 2) 



that can be used to reduce v A to one, at the expense of increasing uq and v 7 . If, on the other 
hand, we eliminate z/ 4 4 + , we obtain the symmetric equation that can reduce z/5 to one. At 
this point, all the powers of the propagators except z/ 2 , z/ 6 and v 7 have been reduced to one. 
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In the same spirit we can derive 

stz/ 2 2+z/ 6 6+ = (2z/ 1567 + v 2 + 2 - 2D) sz/ 6 6+ + (z/ 467 + 2z/ 5 - D) su 2 2+ 

— 2(D — v m7 - 2u 5 ) (z/ 6 6+7- + z/ 4 4 + 5~ + u 3 3 + l~ + ^2+1") 

+ s (2z/ 3 3+ + z/ 2 2+) [(V 7 7+ + ^ 6 6 + + M + ) 5 - - z/ 7 7+l _ " 

+ tz/ 2 2+z/ 6 6+ (3- - I - ) + 2sz/ 4 4+z/ 6 6+5- 

+ 2(D- Z/ 467 - 2z/ 5 ) (2D - 2Z/ 157 " ^2346 ) , (3.47) 

that, together with the symmetric one for i/ 2 2 + is 7 7 + and with 

su 6 6 + u 7 7 + = (D - z/ 567 - 2z/ 4 - 1) z/ 6 6+ + (D — z/ 467 - 2u 5 - 1) u 7 7 + 

+ z/ 6 6+z/ 7 7+ (3- + I - - 5- - 4-) - z/ 4 4+z/ 7 7+5- - u 5 5 + u 6 6 + 4~ 

+ z/ 7 (z/ 7 + 1)7++ (l- - 5-) + v 6 (v 6 + 1)6++ (3- - 4-) , (3.48) 

reduces all powers except one [y 2 or z/ 6 or u 7 ) to unity. 

We can decrease v 2 at the expense of increasing uq and v 7 using 

[(u4 -v 7 + 2z/ 23 + 2- D) s + (i/45 - z/ 67 ) i] z/ 2 2 + = (z/ 5 - 1/3) sz/ 4 4 + + (z/ 7 - z/ 3 ) sz/ 6 6 + 

- (D — 2v 7 - z/ 16 - 2) su 7 7 + — (D — 2v 5 - u u - 2) sz/ 5 5+ 
+ (t + s)u 2 2+ [v 7 7 + (l- - 5-) + z/ 4 4+ (6- - 3-)" 

+ tz/ 2 2+ [z/ 5 5+ (7- - 1") + z/ 6 6+ (3- - 4-)" 

- sv{\ + (z/ 7 7+5 - + z/ 5 5 + 7~) 

+ sz/ 3 3+ [(2z/ 7 7+ + z/ 4 4+) 6 _ + (2z/ 5 5+ + z/ 6 6+) 4" 

- {2D - 2u 57 - 3z/ 46 ) [u 6 6 + 7~ + z/ 4 4+5~ + (z/ 3 3+ + z/ 2 2+) 1 
+ (i/57 - 2i/ 2 ) [z/ 7 7+6- + z/ 5 5+4- + (z/ 2 2+ + ^1+) 3~" 

+ 4L> 2 - 2 (5z/ 57 + 4z/ 46 + 1/3 - z/ 2 + 2i/i) D + v 7 (5v 17 + 10z/ 456 + 4z/ 3 + z/ 2 ) 

+ z/ 6 (3z/ 6 + 10z/ 5 + 6z/ M + 3z/ 3 - z/ 2 ) + u 5 (5z/ 5 + 10z/ 4 + 4z/ 3 + z/ 2 + 5z/ x ) 

+ v 4 (3^4 + 3z/ 3 - i/ 2 + 61/1) - 2v 2 (2i/ 3 + z/ 42 ) . (3.49) 



The power of the seventh propagator can be reduced with 

v , ,, b m w , , ,,(£>- 6) f + (5£>- 21/7- 26) s„. 

iJ — zz/ 7 — o 
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(2D-u 7 -ll)t + (3D-2u 7 -15)s, 

+ p o v 7 ( 1 

D — v 7 — 5 



(D - 2u 7 - 4) t + (5 - D) s A , . „ , , r , 
+ p- — - K - ^4+6+3- + pit + s)\-v 7 {v 7 + 17++6- 

2 (3~ - 6~) + 1 



L> - v 7 - 5 
+ 5+z/ 7 7+ (3- -4- - 6~) + 4+z/ 7 7+ 

- 2+ (z/ 7 7+ + 4+) 1~ + 5+6+3- - 2 5++4 - } 

- 2 ^ n ^ g . (4+z/ 7 7+5- + 5+6+4-) + (t + a) [2 6++ (4" - 3") 

iJ — 2z/ 7 — o v ' L v y 

+ v 7 {y 7 + 1)7++ (5- - I - ) + 6+z/ 7 7+ (5- + 4- - 3- - l - )" 

3+5+ (I - - 4- - 7-) - 3+6+ (4- + 7-) - 3+z/ 7 7+ (4- + 2~) 



+ ps 



4 + ( z/ 7 7+2 _ + 6 + 7 _ + l + 2 _ + 3 + 7 



+ (2D - 3v 7 - 7)p[6+7 _ + 4+5- + (3+ + 2+) 1~ - 2(D - v 7 - 4)] , (3.50) 



where we have introduced the shorthands 
D - 6 



P 



a 



D — 2v 7 — 6 

(5D 2 - 8u 7 D - 50D + 2z/f + 42z/ 7 + 124) s + (2D 2 - 3z/ 7 L> - 21D + 18u 7 + 54) t 



D — 2z/ 7 — 6 

Equation (|3.50|) is not as general as the previous ones since we have set all the powers of the 
other propagators to unity. In addition, since this equation contains 7++, we cannot always 
reduce v 7 to one, but we also have integrals where v 7 — 2. A similar identity can be obtained 
by symmetry for 6++, so that we are left with three integrals: Xbox"°(l, 1, 1, 1, 1, 1, 1; s,t), 
Xbox D (l, 1, 1, 1, 1, 1, 2; s, t) and Xbox D (l, 1, 1, 1, 1, 2, 1; s, t). 

The last step is to write the integral with u e = 2 as a combination of the other two. This 
can be done with the identity that links 6+ with 7+. We derived such an identity, equating 
the expressions obtained by acting with i/ 7 7+ on ^2+^44+ and by acting with z^2+ on 
z/ 4 4+z/ 7 7+ 

(D - 6)(D - 5)— 6+ = (D — 6)(D — 5)7+ - 4 (D ~ 5)3 
t I s t I s 

- (5+7+ + 1+7+ + 1+5+ + l+4+)6- - (3+7+ + 4+6+ + 3+6+ + 3+4+)5~ 

+ (3+4+ + 4+6+ + 5+6+ + 3+6+) I - — (D — 7)4+7+1- 
1 

- -(2+7+6- + 2+5+6- + 4+7+5- + 2+6+5- - 1+4+3" 



t + s 



2+4+7+ (6~ + 5" 



2 (2++7 + - 2 + 4 + 7 



+ -7+ 



}++A + 
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- -{2+ [3+ (6+ + 5+ + 4+) 7- + 4+6+7- + 3+ (7+ + 6+ + 5+) 4" 
+ 3+5+ (2 7+ - 2+) !-}-(£>- 6)|±j^j [2+ (7+4" + 4+7~ + 5+4") 



+ 2(l+7+ + 5+7+ + 1+6+ + l+5+)4 _ - 2 1+6+3" 
- (D - 6) [4+7+3- + (3+4+ + 4+6+ + 3+6+ + 3+5+) 7" 
+ (D — 6)—^— [7+3+4- + (7+3+ + 7+4+ + 1+4+) 2" 



2(D - 5)t + s 



+ (£>-5) 



+ (£>-5) 



(D-5) 



t + s 

(£>-5)f + (2£>-ll)s 



2s 



2D - 13 



2+7+5" 



4+7+6" 



s(t + s) 



(2+ + 1+) 3- + 7+6" 



t + s 



2 3+7+6- + 3+4+6- - 1+7+5 - - 1+5+7" 



-2+5+3" 

2s 



(D - 5)t 



s(t + s) 

+ (D - 5) [I - (3+5+ + 2 2++ + 2+3 



4+5- + 2+1- + 6+7- + 3+1" 



2D - 15 



3+7+1" 



+ 



2+ (5+7+1- + 4+6+3-) 



I 5*-2(£>-7) V 7+1 . 



(D-Q) 



2(t + s) 
t + 2s 



2s 



2+7+3" 



(2D - 9)t 2 — (D — 5)st - 2(D - 5)s s 
2s(i + s) 

2(D - 5)t 2 + (D — 6) st + 2(D - 6)s s 
2s(t + s) 



2+5+1" 



2+6+4" 



+ (SD-W)t + 2(D-5)s 2+7+1 _ + (D-4)t + 2s 2+6+r 



+ 



2s 



2(D - 5)t 2 — (D — 6)st - 2(D - 6)s 2 



2s(t + s) 



2(t + s) 
2+5+7- - (4£>-21) : 



(D-5)t+(D-6)s 3+6+i , 

t + s 



2(f + s) 

QD-4)t-2p-6)s 2+i+3 _ 

2s 



4+6+3" 



+ (2D - ll)t 2 + (D - 7)st + 2{D - 6)s 2 2+6+3 _ + f + 2£ 2+6+r 



+ 



2s(t + s) 
2(D - 5)* + (D - 6)s 



2s 



t + s 



3+5+4- + (D — 



(D-5)t+(2D-ll)s 
^ s(t + s) 



5+4" 



+ ^^2+4+5- - — — — — ^ S 2+4+l~ + 2 ^- 5 ) t - g 2 +4+6" 

2 2s 2s 

where we have set all the powers of the propagators to unity. 



(3.51) 



At the end of this reduction program, we are left with the following two crossed-box 
integrals: Xbox D (l, 1, 1, 1, 1, 1, 1; s, t) and Xbox D (l, 1, 1, 1, 1, 1, 2; s, t), plus simpler diagrams 
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that can always be expressed as a combination of: 

- the master crossed triangle of Fig. ^| (a) 

Xtri D (s) = Xbox D (1, 0, 1, 1, 1, 1, 1; s, t) , (3.52) 

- the master diagonal box of Fig. (b), produced by 

Dbox D 0, t) = Xbox D (0, 1, 1, 0, 1, 1, 1; s, t) = Xbox D (1, 1, 0, 1, 1, 1, 0; s, t) , (3.53) 
together with 

Dbox D (s, u) = Xbox D (1, 1, 0, 1, 0, 1, 1; s, t) = Xbox D (0, 1, 1, 1, 1, 0, 1; s, t) , (3.54) 
and 

Dbox D (t, u) = Xbox D (0, 1, 0, 1, 1, 1, 1; s, t) , (3.55) 

- the master box with a bubble insertion of Fig. || (c), produced by 

Bbox D (s, t) = Xbox D (1, 1, 1, 0, 1, 1, 0; s, t) , (3.56) 

together with 

Bbox D (s, u) = Xbox D (1, 1,1,1, 0, 0, 1; s, t) , (3.57) 

- the master triangle with a bubble insertion of Fig. |2| (d), produced by 

Btri D (s) = Xbox D (1, 0, 1, 0, 1, 1, 0; s, t) = Xbox D (1, 0, 1, 1, 0, 0, 1; s, t) , (3.58) 

- the master sunset diagram of Fig. ^| (e), produced by 

Sset D (s) = Xbox D (0, 0, 1, 0, 1, 0, 1; s, t) = Xbox D (1, 0, 0, 1, 0, 1, 0; s, t) , (3.59) 
together with 

Sset D (t) = Xbox D (0, 1, 0, 0, 1, 1, 0; s, t) , (3.60) 

and 

Sset D (u) = Xbox D (0, 1, 0, 1, 0, 0, 1; s, t) , (3.61) 

where we explicitly used the symmetry of the crossed box 

Xbox D (z/i, zaj, z/ 3 , z/ 4 , z/ 5 , u 6 , v 7 ; s, t) = Xbox D (v 3 , z/ 2 , v x , u 7 , v 6 , u 5 , z/ 4 ; s, t) . (3.62) 

The reduction and dimensional shift of the generic scalar diagonal box and of the box with a 
bubble insertion have been treated in Ref . [[| . The tensor reduction and dimensional shift of 
the crossed triangle is treated in Appendix [A]. All the diagrams obtained from the pinchings 
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(a) (b) 




(e) 

Figure 2: The master integral topologies of the pinchings coming from the reduction 
of the scalar crossed two-loop box. 



of the other propagators can be easily reduced to the previous ones using the integration- 
by-parts identities for the crossed box and for the crossed triangle, with the corresponding 
powers of propagators set to zero. 

Instead of keeping Xbox D (l, 1, 1, 1, 1, 1, 2; s, t) as one of the two members of the basis, we 
prefer to switch to a more symmetric integral: Xbox Z) (l, 2, 1, 1, 1, 1, 1; s, t). The expression 
for this one, in terms of the other two, can be easily obtained through the application of the 
reduction formalism outlined above 

Xbox D (l, 2, 1, 1, 1, 1, 1; s, t) = a Xbox D (l, 1, 1, 1, 1, 1, 1; s, t) 

+ c 2 Xbox D (l, 1, 1, 1, 1, 1, 2; s, t) + c 3 Xtri D (s) + c 4 Dbox D (s, t) 

+ C5 Dhox D (s, u) + cq Dbox D (t, u) + C7 Bbox D (s, t) + c$ Bhox D (s, u) 

+ c 9 Btri D (s) + C10 Sset D (s) + c u Sset D (t) + c 12 Sset D (u) , (3.63) 

where the coefficients Cj are collected in Appendix |B[ This allows us to define as master 
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integrals 



Xboxf (s,f) = Xbox D (l,l,l,l,l,l,l;s,t) 
Xboxf 0,t) = Xbox D (l,2,l,l,l,l,l;s,t), 



(3.64) 
(3.65) 



which are symmetric under the exchange t <-> u. This will produce more compact results in 
the rest of the paper. 



3.2 The dimensional shift for the two master integrals 



Exploiting the equivalence between d and V of Eq. (|3.19|) , we can rewrite Eq. (|3.20 ) 
using (PI) 



Xbox D (pi, z/ 2 , z/ 3 , z/ 4 , z/ 5 , u 6 , u 7 ; s, t) 

= ' (z/ 7 7+ + z/ 6 6+ + u 5 5+ + u 4 4+) (z/ 3 3+ + z/ 2 2+ + ^1+] 

+ (z/ 5 5 + + z/ 4 4+) (z/ 7 7+ + z^ 6 6 + ) Xbox D+2 (z^, u 2 , z/ 3 , z/ 4 , z/ 5 , z/ 6 , u 7 ; s, t) , 

(3.66) 



and more specifically 

Xboxf (s, t) 

Xboxf (s, t) 



(7+ + 6" 



5+ + 4+) (;r • 2' 

D+2 



+ (5+ + 4+) (7+ + 6+) j Xboxf +2 (s, t) 
5+ + 4 + ) (3 + + 2 2+ + 1+ 



(7+ + 6" 



+ 5 



Xboxf +2 (s, *) . 



(3.67) 



(3.68) 



The right-hand side of the two equations can then be reduced to a linear combination of mas- 
ter integrals in D + 2 dimensions, and the system can then be inverted to give Xboxf +2 (s, t) 
and Xboxf" 1-2 (s, £) as a function of master integrals in D dimensions 



where 



Xboxf +2 (s,t) = A(t,u)+A(u,t), 
Xboxf +2 (s,t) = B(t,u) + B(u,t), 

A(t, u) = a x Xboxf (s, t) + a 2 Xboxf (s, i) + a 3 Xtri D (s) + a 4 Dbox D (s, t) 
+ a 5 Dbox D (t, u) + a 6 Bbox D (s, t) + a 7 Btri D (s) 
+ as Sset 15 (s) + ag Sset D (i) , 



(3.69) 
(3.70) 
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B(t,u) = &i Xboxf (s,t) + 6 2 Xboxf (s,t) + 6 3 Xtri D (s) + 6 4 Dbox D (s, t) 
+ b 5 Dbox D (t, u) + & 6 Bbox D (s, t) + b 7 Btri D (s) 

+ 6 8 Sset D (s) + 6 9 Sset D (t) . (3.71) 
Some of the coefficients a« and 6j of the system are collected in Appendix y. 

4 Differential equations for the two master integrals 

In the previous sections we showed how tensor integrals can be expressed in terms of the 
two master crossed boxes, Xboxf and Xboxf plus simpler master diagrams. The analytic 
expansion in e = (4 — D)/2 of the first master integral was computed in Ref. 0. We can 
obtain the analytic form for the second one by writing the derivative of Xboxf with respect 
to one of the two independent physical scales (that we choose to be t), as a combination of 
master integrals, and solving the equation for Xboxf. 

Moreover we can verify the correctness of both the expressions of Xboxf (s, t) and 
Xboxf (s, i), by deriving an analogous differential equation for Xbox^*, and checking that 
the obtained identity is satisfied. 

Using the procedure outlined in Section ^ we differentiate with respect to t the two 
master integrals written in the form of Eq. (|3.12p . In this way, the only i-dependence comes 
from Q of Eq. ( p. 15 ) and the result is 



d 

— Xboxf (s, t) = Xbox D+2 (1, 2, 1, 2, 1, 1, 2; s, t) - Xbox D+2 (1, 2, 1, 1, 2, 2, 1; s, t) , 
d 

— Xbox^ (s, t) = 2 Xbox D+2 (1, 3, 1,2,1, 1, 2; s,t)-2 Xbox D+2 (1, 3, 1, 1, 2, 2, 1; s, t) 



Applying the reduction formulae of Section |3.1| and the dimensional-shift of Section |3.2j to 



the right-hand sides of the system, we can rewrite them as a combination of master crossed 
boxes and pinchings in D dimensions. 

The final system of differential equations in arbitrary D is given by 

|-Xbox?(M) = -^—[H(t,u)+H(u,t),} (4.1) 
at t — u 

|-Xbox£( S ,t) = J—[K(t,u) + K{u,t),) (4.2) 
at t — u 

where 

H(t,u) = hi Xboxf (s, t) + h 2 Xboxf (s, t) + h 3 Xtri D (s) + h 4 Dbox D (s, t) 
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+ h 5 Dbox D (t, u) + h 6 Bbox D (s, t) + h 7 Btri D (s) 

+ h 8 Sset D (s) + hg Sset D (t) , (4.3) 
K(t, u) = ki Xboxf (s, t) + k 2 Xboxf (s, t) + k 3 Xtri D (s) + k 4 Dbox D (s, i) 
+ k 5 Dbox D (t, u) + k 6 Bbox D (s, t) + k 7 Btri D (s) 

+ k 8 Sset D (s) + k 9 Sset D (t) , (4.4) 
and the coefficients hi and ki are listed in Appendix [D]. 

From the symmetry t <-> w of the two master integrals, we expect the two derivatives to 
be anti-symmetric with respect to the exchange t <-> w, and this is a further check of the 
correctness of the system. 



5 The off-shell method 



An alternative way to derive the system of equations (|4.1|) and (|4.2|) is based on the con- 
struction of a set of differential equations for the crossed-box integrals with the momentum 
Pi taken off-shell. This is done using the algorithms and computer programs described in II , 
by means of which the differential equations for any massless two-loop four-point function 
with three light-like and one off-shell leg can be obtained. We then take the on-shell limit 

pI^o. 

To derive the differential equations for the diagrams with one leg off-shell, we take p 2 , Pz 
and p4, all massless and on-shell, as independent momenta, and form the derivatives with 
respect to the external invariants S23, S24, S34, where —p\ 

9 1 f, , 9 , , 9 M 9 \ 



'ds 23 2 V Opt} ™ dp% ™ dp% 

9 1 ( , n 9 v 9 I 9 \ ( ~ ~\ ° 

824 dsV, = 2\ + ^M~ p 'M +fi M^ 



9 l/^d d ^x 



The derivatives with respect to p% act on the propagators in the Feynman integrand of 
Eq. ( |2.ip and give rise to a number of different integrals with increased powers of the pro- 
pagators. Using integration-by-parts and Lorentz-invariance identities, all these different 
integrals can be reduced to a small number of master integrals. In particular, any crossed- 
box topology integral can be reduced to a linear combination of two master crossed-box 
integrals, which we can choose to be the off-shell versions of Xboxf and Xboxf, plus a num- 
ber of pinched master integrals. Thus, we obtain a system of two coupled linear first-order 
differential equations for the two master crossed-box integrals. 



17 



Since, in the on-shell limit, we are interested in the derivatives at fixed pf = 0, we rewrite 
the system of equations in terms of the variables pi = S23 + S24 + S34, s = S34, t = S23, and 
the corresponding derivative operators: 



a 


d 


d 


ds 


<9s 34 


8s 2 4 





d 


d 


Of 


ds 23 


8s 2 4 


d 


d 




dpi 


ds 2 4 





(5.2) 



The next step is to take the on-shell limit p\ — >0. A complication here is the fact that factors 
of 1/pl and l/{pl) 2 appear in some of the coefficients of the pinched master integrals. It is 
also important to realize that several of the off-shell master integrals with less than seven 
propagators become reducible in the on-shell limit (an example of such a case will be given 
below). When these reductions are taken into account, all the terms proportional to 1/pl 
and l/(pl) 2 cancel out, so that the on-shell limit of the system of equations for the two 
master crossed-box integrals is indeed well defined. 

However, because of these factors of pi in the denominator, it is not sufficient merely to 
replace all pinched master integrals by their limit pi — >0: subleading terms in pi also have to 
be included. The Taylor series around pi = for the pinched diagrams can be obtained in 
a straightforward manner by using the differential equation in pi. As an example, we con- 
sider the propagator diagram Sset D (pi — s — t), which fulfills the homogeneous differential 
equation 

4* Sset D (pi - s - t) = 9 D ~ 3 Sset D (p\ — s — t) . (5.3) 
dpi v ' V\ — s — t v ' 

Iterating this equation, we obtain 

Sset D (^_ s _t) = Sset D (-s-t)+pl^-^- t Sset D (-s-t) 

+ \v\ {D ^ (D t) ~2 A) ** D (- - *) + O (p?) • (5.4) 

If a factor of 1/pl appears in the homogeneous term of the differential equations for an 
off-shell master integral, this integral becomes reducible in the limit p\— >0. An example 
of a reducible master integral is the vertex diagram Rtri(pf,s) of Fig. |3|, which fulfils the 
differential equation 

d r, ■ ( 2 \ D ~ 4 2 Pl - s r, • / 2 \ 3D - 8 1 n D / \ 

dpj Rtn ^ s ) = — pM^V) Rtn s ) - pM^s) {s) ' ( } 

Multiplying by p\, one can immediately read off the limit 

Rtri (p?-0, s) = \ Sset D (s) + O (pi) . (5.6) 
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Figure 3: The master integral Rtri(p^,s)j which becomes reducible when p\ is 
on-shell. 



Subleading terms in p\ are again obtained by successive iteration of the differential equa- 



All the information required to derive the massless limit of the system of equations for 
the master crossed-box diagrams is contained in the differential equations for the pinched 
master integrals. 

The differential equations obtained with this method are in perfect agreement with ( |4.1|) 
and ( |4.2|) , giving an independent confirmation of the correctness of the two methods. 

In addition, we checked that the two methods give the same results for the crossed boxes 
with the powers of propagators equal to unity and with one scalar product in the numerator, 
of the type pi ■ k or pi ■ I. This was done computing Eq. (|3.21|) with the method described in 
Section |3| and then contracting the result with one of the external momenta. 

6 Analytic expansion of the second master integral 

Inserting the e expansion of Xboxf (s, t) computed in Ref. || and the e expansions of 
the sub-topologies listed in Refs. || £|, [□], into Eq. ( |4.1| ), and solving it with respect to 
Xbox® (s,t), we obtain, in the physical region s > 0, t,u < 0, 



tion QOP . 



Xbox^(s,t) =r 2 (l + e) 



( 



Gi(t,u) G 2 (t,u) G x (u,t) G 2 (u,t) 



s 3 t s 2 t 2 S 3 U S 2 V? 



(6.1) 



where 
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1 8 
+ - f 1 - 12 7r 2 - 24 T + T 2 - 24 U + 16 T U + U 2 ) - 43 - 18 T + 13 T 2 + - T 3 
e v / 3 

8 / 112 

- 18 U + 16 T C/ + 11 T 2 U + 13 P 2 - 20 T f/ 2 + - U 3 + ix 2 17 T + 17 C/ 

3 V 3 

- 122 C(3) + 62 T Li 2 f--^ - 62 Li 3 + 62 S 1)2 (~ 



+ Z7T 



(16 + 6 T + 6 C/) - 34 - 9 n 2 - 6 T - 10 T 2 - 6 E/ + 14 T J7 - 10 U 2 



(6.2) 



-2c 



G 2 (t,u) = s 

1 / 29 5 



e 4 e 3 



2 2 



+ 



+ 



12 



n 2 + 7T _ T 2 + 2 QU -4TU -U 2 



T 3 



7T 



— + 17T + 2T 2 + — (14 + 5 T - 29 £/) + 13 £/ - 28 T P - 4 P 2 

2 3 6 V ; 



2 ' ' 19 = ^ ^ • 2 Li ' ' 



+ 3TP 2 -f/ 3 + — C(3) -2TLi 2 

2 v s 



+ E + E n * + 7T -5T 2 ~—T 3 + -T A + 5U -2QTU + -T 3 U -2U 2 
2 40 3 3 3 

4 4 

+ 2AT U 2 — T 2 U 2 — 8 U 3 TU 3 + - U 

3 3 



+ — (79 - 22 T - 5 T 2 - 200 C/ + 76 T £/ + 25 U 2 ) + (68 - 13 T - 33 [/) C(3) 
5 v / 

+ (10tt 2 -32T + 17T 2 + 12 TU) Li 2 + (32 - 60T - 12 U) Li 3 / / 

+ (28 T - 6 U - 32) Si, 2 f-^) - 26 S lj3 - 36 S 2 , 2 (--^ + 86 Li 4 



+ m 



1 



31 



— + — (11 - T + U) + - 1 7r 2 -10T-2T 2 + 4[/-2T[/-2^ 

e 3 e J e V 6 



10 



7T 



+ 11+4T-2T 2 + — T 3 + — (-65 + 28T-P) + 2P-8TP-8P 2 



-2P :! - S!)U3) + (I J V+ 1ST) Li, (-^ -32Li 3 (-^)+44S 1)2 (-- 



(6.3) 



and T = log(— t/s), U = log(— u/s). We used Nielsen's generalized polylogarithms S niP 
14fl , defined by 



Sn,p (-^) 



^n+p-l i l 0g »-l( t ) l 0g P(l _ ^ 



(n — l)\p\ Jo 



dt 



n,p > 1, x < 1, 



(6.4) 
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Figure 4: The physical regions (i), (ii) and (iii) in the (s,t,u) -plane. 

where the usual polylogarithms are given by 

Li n (x) = S„_i,i (x) . (6.5) 

The three kinematically accessible regions of the phase-space are depicted in Fig. 0L 



(i) s > 0, t, u < 0. All logarithms and polylogarithms occurring in Eqs. (|6.2j ) and ( 3.3 
are real in this region. 

Formulae for the other two regions, (ii) and (iii), can be derived by analytic continua- 
tion, starting from region (i) and following the paths indicated in the figure. 

The analytic continuation can be performed through a few simple steps. 

(ii) t > 0, s, u < 0. Going from region (i) to region (ii), we have to pass through two 
branches: t = and s = 0. We can then split the analytic continuation into two steps: 

- we first split the logarithm T = log(— t) — log(s). At t = 0, nothing happens to the 
polylogarithms S n>J , (—t/s), but log(— t) gets an imaginary part: log(— t)— > log(t) — 

Z7T. 

We are now in an unphysical region, where both s and t are positive and u is 



negative. Using the transformation formulae for x— >l/x (see, eg. Refs. [ 14], [3J), 
we can express S„ iP (—t/s) in terms of S„ )P (s/t) and log(t/s). 

To enter region (ii), we have to pass now the branch point at s — 0. We split 
log(t/s) = log(t) — log(s) and U = log(s + t) — log(s) and we analytically continue 
log(s)^log(-s) + m. 
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In this way, for example, the logarithms in Eqs. ( |6.2| ) and ( |6.3| ) undergo the transfor- 
mation 

log(-£) bg(~) -2i7r, (6.6) 

l0g (~s) ~* log (s)~ i7r ' ( ' 6 ' 7 ' ) 

(iii) w > 0, s, t < 0. The procedure to go from region (i) to region (iii) is similar to the 
previous one, but it requires an additional step. 

- We rewrite S n>p (— t/s) in terms of S„ )P ((s + t)/s), log(— t/s) and log((s + t)/s), 
using the transformation x — > 1 — x, and we split the logarithms as before. In 
passing the first branch point at u — 0, the poly logarithms are well defined while 
log(s + t)—> log(— s — t) — in. 

- We invert now the argument of the polylogarithms, expressing S n p ((s + t)/s) in 
terms of S niP (s/(s + 1)) and log((— s — t)/s) = log(— s — t) — log(s). Finally, 
log(s)— > log(— s) + in, as we pass the branch point at s = and enter region (iii). 

The logarithms in Eqs. ( p.2| ) and ( |6.3| ) undergo the transformation 

log(-£) - logQ)--, 

l0gf-^ -> log^-^-2z7T. 

The expression for G\(t,u) and G2{t,u) in this region can also be obtained directly 
from the expressions in region (ii), using the symmetry t <-> u. 

A non-trivial check of the correctness of the expressions of Xboxf (s, t) and Xbox^ 3 (s, i) 
comes from Eq. ( f4.2| ), that must be identically satisfied, once the respective e expansions are 
used. 

7 Conclusions 

High-energy scattering processes are one of the most important sources of information 
on short- distance physics. Recent improvements of experimental measurements demand 
nowadays the knowledge of 2 — > 2 scattering rates at two-loop order. 

Two non-trivial topologies characterize these processes: planar two- loop boxes 0, |3[], and 
crossed two-loop boxes. 
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(6.8) 
(6.9) 



The crossed boxes constituted the last barrier towards the completion of this goal. In this 
paper we presented an algorithm to reduce tensor two-loop massless crossed boxes with light- 
like external legs to two master crossed boxes (the first one with all powers of propagators 
equal to one, and the second one with the power of the second propagator of Fig. H equal to 
two), plus simpler diagrams. 

We derived the equations that connect the two master integrals in dimensions D and 
D + 2 and the system of first-order coupled differential equations satisfied by the two master 
integrals. 

This last part was done following two different methods: using the raising and lowering 
operators, and using the on-shell limit of the differential equations for the crossed box with 
one off-shell leg |§. The agreement between the results of the two methods is a strong 
support for the validity of the two different procedures. 

The differential equation for the first master integral allowed us to derive an e expansion 
for the second master integral, once the known || expression of the first master integral 
is inserted. The differential equation for the second master integral then provided us a 
non-trivial way to check the obtained expansion. 
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A Tensor reduction and dimensional shift of the crossed two-loop 
triangle 

For completeness, we give here the reduction formulae for the crossed two-loop triangle 
that appears as a sub-topology in the scalar reduction of the crossed box. We make use of 
integration-by-parts and Lorentz-invariance identities to reduce the generic scalar triangle to 



one master integral which has all powers of propagators equal to one ITT], [l2[ , plus pinchings. 



For an alternative solution to this reduction problem, exploiting a connection with massless 
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Figure 5: The generic two-loop crossed triangle. 



three-loop propagator integrals, see Ref. |L5 . 

The generic two-loop scalar crossed triangle of Fig. |5] is given by 

d D k r d D l 



Xtri" (i*, u 2 , u 3 , u,, u 5 , ^s)^J — J— A?A?A?ATA?A? , (A.l) 
where the propagators are 

A x = {k + l + p 3 + Pi } 2 +i0, (A.2) 
A 2 = {k + l) 2 + t0, 



A 3 = l 2 + i0 



2 



Ai = (l+psY + iO, 
A 5 = k 2 + iO, 
A 6 = (k + Pi ) 2 + iO. 

The external momenta are in-going, two of them being light-like, p 2 = p\ = 0, while p\ 2 = s 
is the only physical scale. 



Repeating the reasoning of Section |3J], we can build eight integration-by-parts identities 
and a single Lorentz-invariance identity. They depend upon one irreducible scalar product 
in the numerator, that we choose to be (/ ■ p^): 

svil+ + {2D - 2z/ 235 - v im ) - z/ 4 4 + 3 - - ^1+2' - z/ 6 6+5 _ = (A.3) 
sz/ 2 2+ + {2D - 2z/ 146 - z/ 235 ) - z/ 3 3+4- - v 2 2+l~ - v 5 5+6~ = (A.4) 
2 (/ • p 4 ) uxl+ — {D — z/ 24 - 2vi) + z/il+ (2- + 4- - 5-) 

+ v 2 2+ (3- - 5 _ ) + z/ 4 4 + 3 - = (A.5) 

2 {I • p 4 ) u 2 2+ — (D — z/i 2 5 - 2u Q ) + v 2 2+ (6~ - 3 
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4+ (I" - 6-) 



- v x l+ (q~ - 4~) - z/ 5 5+6 - = 

2 (Z • p 4 ) ^3+ + (D - z/345 - f 6 ) + u 3 3 + (2" - 5~) + 1/4 

- z/ 5 5+6- = 

2 (Z • p 4 ) ^4+ + su 4 4 + — (D — z/ 346 - 2u 5 ) + z/ 3 3+ (5 - - 2") 

+ z/ 4 4+ (6 _ - l - ) + z/ 6 6+5 _ = 
2 (Z • p 4 ) z/ 5 5+ + sz/ 5 5+ — (D — u m - 2z/ 4 ) + z/ 5 5+ (4- + 6~ - l - ) 

+ z/ 6 6+ (4- - 1") + z/ 3 3+4- = 
2 (Z • p 4 ) z/ 6 6+ — {D — z/ 45 - 2z/ 3 ) - z/ 6 6+ (3" + 5" - 2") 



(A.6) 



(A.7) 



(A.8) 



(A.9) 



+ z/ 5 5+ (2- - 3 _ ) + z/ 4 4+3 - = (A.10) 

2 (Z • p 4 ) - (£> - z/ 2356 ) + (2- + 6- - 5-) + z/ 4 4+3- = 0. (A.ll) 

By eliminating the irreducible scalar product in the numerator, we obtain 

= -{2D- 2z/ 235 - u 14s ) + z/ 4 4 + 3 _ + i/il+2 _ + z/ 6 6 + 5 _ , (A.12) 

which, together with the symmetric one for z/ 2 2 + , can reduce v\ and 1/2 to unity. Using 

1 



z/ 3 3 



+ 



D — 2 — z/ 3456 
1 



+ 



D -2- 2z/ 34 



(z/ 4 4 + z/ 6 6+l _ - z/ 3 3 + z/ 5 5+2 

(D — 2 — 2v m ) z/ 6 6+ + 2 (1/3 - i/ 6 ) z/ 5 5+ 



(A.13) 



and the symmetric one for z/ 4 4 + , we can reduce z/ 3 and 1/4 to one. To complete the reduction, 
we use 

(^56 - M = ^2+ (3- - 5-) + (4- - 6-) , (A.14) 

that can be re-iterated until (zv 56 — z/ 34 ) = 0. Since we are applying this identity to scalar 
integrals where u 3 and 1/4 have already been reduced to one, the reduction procedure will 
stop when u 5 — u 6 — 1. This integral cannot be reduced any further, and we choose the 
crossed master triangle to be 



Xtri D (s) = Xtri D (1, 1, 1, 1, 1, 1; s) . 
The application of the identity 



(A.15) 



d~ = (z/ 6 6+ + z/ 5 5+ + z/ 4 4+ + z/ 3 3+) (z/ 2 2+ + ^1+) + (z/ 4 4+ + i/ 3 3+) (z/ 6 6+ + v 5 5 + ) 

(A. 16) 
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to Xtri £)+2 (s) and its further reduction, gives rise to the dimensional-shift formula 

Xtri D + 2 (s) = ^~ A K, rXtri D ( S ) 

v ' 4(D - 2)(2D - 7) (2D - 5) V ' 

37D 3 - 313D 2 + 858D - 752 „ n , , 

Btn (s) 



2(D - 4)(D - 2)(2D - 7) (2D - 5) (3D - 8) 

43D 4 - 478D 3 + 1963D 2 - 3530D + 2352 D 
+ 2(D - 4) 2 (D - 3)(D - 2)(2D - 7)(2D - 5)s {s) ' ( } 



The expression of the master integral of Eq. (|A.15|) has been computed in Refs. |TTL [12 



B Coefficients c, 



We collect here the coefficients q of Eq. ( j3.63| ) . 
A (D - hfu 



c 2 



c 3 



C6 



C 7 



C8 



Cg 



(D - 6)st 

£ — u 
t 

(P-4)(2P-9) 
(D - 6)st 

(D — 4) (3D - 14)w 2 
~ 3 (D-6) 2 (2D-ll)s 3 i 3 

(D-4)(3D-14)t 
(D-6) 2 (2D-ll)s% 2 

(D - 5)(D - 4) (3D - U)s [2(D - 5)t - (D - 6)s 
(D-6) 2 (2D-ll)t% 2 

(D - 3) (3D - 14) [(5D - 26)t + 2(D - 5)s] 
(D - 6) 2 st 3 

(D - 3)(3D - 14) [(5D - 26)t + 2(D - 5)s] 



2 (7D 2 - 76D + 206) t + (D — 5)(5D - 28)s 
2 (7D 2 - 76D + 206) t + (D - 5)(5D - 28)s 



(D - 6) 2 s^ 2 
(D - 3) (3D - 10) 



(41D 3 - 620D 2 + 3124D - 5248) st 3 



ClO 



2(D-6) 2 (D-5)(D-4)s 3 t% 2 
+ 2(D - 6) 2 (D - 4)t 4 - (61D 3 - 922D 2 + 4640D - 7776) sH 

- 4(D - 5)(3D - 14)(4D - 21)s 3 t - 4(D - 5) 2 (3D - 14)s 4 

(D-3)(3D-10)(3D-8) 
~ (D - 6) 2 (D - 5)(D - 4) 2 (2D - ll)s 4 t 3 w 2 
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Cll 



Cl2 



x [2(D - 6)(D - 4) (2D - 11)(5D - 22)t 4 

- 2 (64D 4 - 1305D 3 + 9978D 2 - 33902D + 43180) st 3 

- (208D 4 - 4209D 3 + 31899D 2 - 107314D + 135216) s 2 t 2 

- (D — 5) (3D - 14) (34D 2 - 359D + 944) s 3 t - (D - 5) 2 (3D - 14) (5D - 28)s 4 

(D - 3) (3D - 14) (3D - 10) (3D - 8) 
(D - 6) 2 (D - 5)(D - 4) 2 (2D - ll)sH 4 u 2 

x 2(D - 5) (7D 2 - 76D + 206) t A + (61D 3 - 952D 2 + 4941D - 8528) st 3 
+ (92D 3 - 1417D 2 + 7253D - 12336) s 2 t 2 

+ (56D 3 - 853D 2 + 4317D - 7258) s 3 t + 2(D - 5)(2D - 11)(2D - 9)/ 

(D — 3)(3D — 14)(3D — 10)(3D — 8) r , , , 9 n 

2(D - 5) (7D 2 - 76D + 206) t 



(D - 6) 2 (D - 5)(D - 4) 2 (2D - ll)s 3 t% 3 
(9D 3 - 130D 2 + 615D - 948) st 3 - 3(D - 6)(D - 5)(D - A)s 2 t 2 

(D - 5) (3D 2 - 30D + 74) s 3 t + (D — 6)(D — 5)V 



C Coefficients of the dimensional-shift system 



:b.i: 



We give here only the first three coefficients of the system ( |3.71|) . The whole list can be 
obtained from the authors (C.A. and CO.). 



«2 



a 3 



bi 



3 (D - 4)s 2 [(3D - 14) (3D - 10)s 2 - 4(5D 2 - 39D + 74)tu] 
512 (D-3) 3 (2D-9)(2D-7)tu 

(D - 6)(D - 4)s [4(D - 3)tu - 3(3D - 10)s 2 ] 
512(D - 5)(D - 3) 3 (2D - 9)(2D - 7) 

(D - 4)s [4(D - 6)(D - 3)tu + 3(3D - 14) (3D - 10)g 2 ] 
512(D - 5)(D - 3) 3 (2D - 7)tu 

s [3(D - 4)(3D - 14)s 2 - 4(D - 6)(D - b)tu] 
128(D - 3) (2D - 9) (2D - fjtu 

(D - 6) [4(D - 3)tu + 3(D - 4)s 2 ] 



128(D - 5)(D - 3) (2D - 9)(2D - 7) 

3(D - 4) (3D - 14)s 2 - 4(D - 6)(D - 3)tu 
128(D - 5)(D - 3)(2D - 7)tu 

27 



(C.l) 



D Coefficients of the differential equations for the two master 
integrals 



The coefficients of the system of differential equations ( |4.1|) and ( f4.2|) for the two master 
integrals are given by 

(D -A)s 2 - Atu 



h - 
h 2 = 

h 3 - 

h± - 

h 5 -- 

h -- 

h 7 -- 

h 8 -- 

hg -- 

and by 

h 
k 2 

h 



Atu 
(D-6)s 
~A(D-5) 

(D — 4) (2D — 9)s 
A(D - 5)tu 

3 (D -4) (3D - lA)u 

2 (D - 5) st 2 

3 (D — 4) (3D - 14)s 2 

2 (D - Q)t 2 u 2 

(D-3)(3D-14) 
(D — 5)t 2 

3 (D - 3) (3D - 10) [(3D - 14) (u 2 + t 2 ) + 2(D - A)tu] 
A (D - 5)(D - A)st 2 u 2 

3 (D - 3) (3D - 10) (3D - 8) [(D - 5) (3D - 14) (u 2 + t 2 ) - (D - 6)(D - A)tu] 
A (D - 5) 2 (D - A) 2 s 2 t 2 u 2 

3 (D - 3) (3D - 14) (3D - 10) (3D - 8 



2 (D -6)(D -5) 2 (D -A) 2 st 3 u 2 
+ (7D 2 - 68D + 164) tu + 2(D - 5) 2 t 2 



(2D - 9)(3D - 16K 



(D.l) 



(D-5) 2 g 
tu 

(D-Q)(u 2 + t 2 ) 
2tu 

(D-4)(2D-9) 
tu 

(D - 4)(3D - 14)w [(5 - D)u + (2D - ll)t] 
6 (D - 6)s 2 t 3 
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k 5 
k 6 
k 7 



(5D - 28)tu + (D — Q)t 2 - 2(D - 5)u 2 
2(D - 6)(3D - U)tu (u 2 + t 2 ) 



ks = 3 



3 (D - 5)(D - 4) (3D - 14)s 3 

2 (D - Q)t 3 u 3 

(£>-3)(3£>-14 ) r 
(£> - 6)st% 

3 (D — 3)(3D — 10) 
2(.D-6)(Z>-4)sW 

- (£> - 5) (3D - 14) (V + t 4 ) + 2 (5£> 2 - 49L> + 118) t 

(g _ 3 )(3D - 10)(3D - 8) r 
(7J-6)(/J-5)(/J-4) 2 S W i 3(D " 5) (3D " 14)tM ( U + * ) 

- (D - 5) 2 (3D - 14) (V + t A ) — (D — A) {jD 2 - 70D + 176) t 

(D - 3)(3D - U)(3D - 10)(3D - 8 



2 u 2 



2 2 
U 



{D-h)\D-2)u i 



(D — Q) 2 (D — 5)(D — 4) 2 s 2 t 4 -u 3 
+ (D - 6) (l3L> 2 - 129£> + 318) tw 3 + 2 (5£> 3 - 80L> 2 + 422£> - 734) t 



2 2 
It 



+ (£> - 6)(D - 5)(5D - 24)t 3 w + (£> - 6)(D - 5) 2 t 4 



(D.2) 
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